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ABSTRACT 

Non-thermal X-ray emission in some supernova remnants originates from synchrotron radiation of uhra- 
relativistic particles in turbulent magnetic fields. We address the effect of a random magnetic field on syn- 
chrotron emission images and spectra. A random magnetic field is simulated to construct synchrotron emission 
maps of a source with a steady distribution of ultra-relativistic electrons. Non-steady localized structures (dots, 
clumps and filaments), in which the magnetic field reaches exceptionally high values, typically arise in the ran- 
dom field sample. These magnetic field concentrations dominate the synchrotron emission (integrated along 
the line of sight) from the highest energy electrons in the cut-off regime of the distribution, resulting in an 
evolving, intermittent, clumpy appearance. The simulated structures resemble those observed in X-ray images 
of some young supernova remnants. The lifetime of X-ray clumps can be short enough to be consistent with 
that observed even in the case of a steady particle distribution. The efficiency of synchrotron radiation from the 
cut-off regime in the electron spectrum is strongly enhanced in a turbulent field compared to emission from a 
uniform field of the same magnitude. 

Subject headings: radiation mechanisms: non-thermal — X-rays: ISM — (ISM:) supernova remnants 



1. EVITRODUCTION 

The forward shocks of supernova remnants (SNRs) are 
the most probable accelerators of cosmic rays to high ener- 
gies. Supporting this idea is the believe that the diffusive 
shock acceleration (DSA) mechanism is capable of acceler- 
ating particles to above 100 TeV in young and middle age 
SNRs. Electrons (and/or positrons) accelerated to ~ 10 TeV 
will efficiently radiate in X-rays if magnetic fields exceed 100 
' fiG. With Chandra 's imaging capability it is possible to dis- 
tinguish the synchrotron structures in the X-ray images of 
• Cas A, RCW 86, Kep l er, Tycho and som e other SNRs (e.g. 
! I Vink & Laming||2003l: iBamba et al.ll2006t iPatnaude & FesenI 
2008h . The non-thermal emission seen by Chandra concen- 
trates in very thin filaments in all SNRs. Bamba et al. (2006) 
also reported on the time evolution of the scale length of the 
filaments and roll-off frequency of the SNR spectra. 

The synchrotron filamentary structures could be due to a 
narrow spatial extension of the highest energy electron pop- 
ulation in the shock acceleration region limited by efficient 
synchrotron energy losses. In this case, strong magnetic 
field amplification in the shock vicinity is required to match 
the observed thickness of the non-thermal filaments (e.g . , 
IVink &Lamingll200l iBamba etal.ll200l IVolk et al.lllool . 
thus supporting the case for highly efficient DSA. An alterna- 
tive interpretation, that the observed narrow filaments are lim- 
ited by magnetic field damping and not b y the energy losse s 
of the radiating electrons, was proposed by lPohl et al.l (l2005h . 

There is an emerging class of SNRs (SN1006, RXJ17 13.72- 
3946, Vela Jr, etc) dominated by non-thermal X-ray emission. 
These SNRs are of particular interest since their X-ray emis- 
sion most likely originates from synchrotron radiation of elec- 
trons accelerated to well above TeV energies. Many of the 
non-thermal X-ray SNRs are known to be radio-faint com- 
pared with most SNRs, and some of these remnants were first 
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discovered in X-rays. The sources mentioned are now known 
to emit TeV photons as well and the origin of this emission is 
of fundamental importance for the origin of cosmic rays. Is it 
inverse-Compton emission from the X-ray emitting electrons, 
or pion-decay emission from proton-proton interactions, or 
both? One way to get the answer is to study the magnetic 
field struct ure in the SNRs. 

Recently lUchiyama et al.l (l2007h reported variability of X- 
ray hot spots in the shell of SNR RXJ1713.72-3946on ~ one- 
year timescales. The authors suggested that this X-ray vari- 
ability results from radiation losses and, therefore, reveals the 
ongoing shock-acceleration of electrons in real time. In order 
for radiation losses to be this rapid, the underlying magnetic 
field would need to be amplified by a f actor of more than 100. 
From multi-wavelength data analysis iButt et all (l2008h con- 
cl uded however, that the mG-scale magnetic fields estimated 
bv lUchivama et al] (l2007h cannot span the whole non-thermal 
SNR shell, and if small regions of enhanced magnetic field do 
exist in RX J1713-3946, it is likely that they are embedded 
in a much weaker extended field. Therefore, it is essential to 
study the effects of magnetic field structure and particle dis- 
tributions on the observed synchrotron emission maps. 

Turbulent magnetic fields, with energy densities approach- 
ing a substantial fraction of the shock ram pressure, are a 
generic ingredient in the efficient DSA mechanism. How- 
ever, the effect of turbulence on synchrotron emission images 
and spectra has not yet been addressed. In this paper, we 
model specific features in synchrotron images arising from 
the stochastic nature of magnetic fields. In ^we describe our 
random magnetic field simulations, and these are used in ^ 
to construct synchrotron emission maps for different frequen- 
cies. Our results show that the rapid time variations seen in 
some SNRs may well be produced by steady electron distri- 
butions in turbulent fields. 

2. SIMULATION OF A RANDOM MAGNETIC FIELD 

In order to simulate numerically an isotropic and spatially 
homogeneous turbulent field we sum over a large number 
(A';,,) of plane waves with wave vector, polarization, and 
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Fig. 1 . — Panels (a) and (b) are simulated maps of synchrotron emission in 
a random magnetic field. Intensities, ■ I(R±,t, u), are shown with a linear 
grayscale (white is maximum). The map (a) is made at 5 keV, while (b) is 
made at 50 keV. The shock position, where the electron distribution peaks, 
is indicated with a dashed line. The stochastic magnetic field sample has 
(B-) = 10""* G and 5 = 1 .0. A dim region "A" is marked to contrast a bright 
region "B" and a dot marked as "Dl". 

phase chosen randomly (see e.g. [ Giacalone & Jokipiilll999l : 
ICasse et alJl200llCandia & Roulelll2004 . i.e.. 

N„, 2 

B(r>0 = 5mA*")(^„)cos(k„-r-c^„(k„) + , (1) 

n=l a=l 

where the two orthogonal polarizations A^°'\k„) (a = 1,2) are 
in the plane perpendicular to the wave vector direction (i.e., 
A'-°'\k„) 1 k„, so as to ensure that V • B = 0). We divided 
k-space between A;,nin = Itt/L^^^ and k,^^^ = In/L^in into 
spherical shells distributed uniformly on a logarithmic scale 
(a number of points in the s-shell is M,). Here L^i„ and L^ax 
are the minimum and maximum scales of turbulence respec- 
tively. The spectral energy density of the magnetic field fluc- 
tuat ions is taken as W(k) oc k~^ , where S is the spectral index 
(c.f. lGiacalone & .Tokipiill999h . The average square magnetic 
field {B^) = J dkW{k) is an input parameter for our model. 
In the particular simulation runs presented below we assume 
(jj„(k„) = Vph • k„. The parameter Vph would be a phase veloc- 
ity in the case of a superposition of linear modes (e.g., mag- 
netosonic waves). In this paper, however, we are modeling 
strong magnetic field fluctuations, so the random Fourier har- 
monics in Eq.([T]i comprising the random field should not be 
exactly associated with any MHD eigen modes. 

Following the prescription given above in Eq.([T]i, we sim- 
ulated the temporal evolution of random magnetic fields fill- 
ing a cube of scale size D, where Lmin = 2 x 1Q~^ttD, Lmax = 
0.2 ttD, N,„ = 4800, and N = 6,Ms = 800. The choice of the 
mode number Nm = 4800 provided the simulated fluctuation 
spectrum to be fitted with power law of index 5=1.0 between 
fcmin and fcmax with accuracy better than 2%. The spectral in- 
dex 6=1 corresponds to the so-called Bohm diffusion model 
in DSA. With the simulated random magnetic fields, we cal- 
culated the local spectral emissivity of polarized synchrotron 
emission at various times and for different distributions of 
electrons or positrons. The technique is described next along 
with the simulated emission maps. 

3. SYNCHROTRON IMAGE SIMULATIONS 

The formation length I / of the synchrotron radiation pattern 
of an individual ultra-relativistic electron is energy indepen- 
dent, i.e., // w Rg/^ = mc^/eB « 1.8 x 10^ B"'^ cm, where 
Rg is the gyro-radius of a particle with Lorentz factor 7. The 
length If is a small (7"') fraction of the particle gyro-radius. 
In this study we only consider the effects of magnetic fluc- 
tuations having scales comparable to the gyro-radii of ultra- 
relativistic particles that are much larger than / /. Therefore, 
the fluctuation wavenumbers k satisfy k-lj <^ 1 . This restric- 
tion allows us to apply the standard formulae for the syn- 
chrotron power of a single particle of Lorentz factor 7 3> 1 



in a homogeneous magnetic field, and then to integrate this 
power over the line of sight through the system filled with 
random field fluctuations. 

We start with the spectral flux densities p\}\9,E) 
and p\^\9,E) with two principal directions of polariza- 
tion radiated by a p a rticle with energy E, given by 
iGinzburg & Syrovatskiil (1 19651) [their Eqs.(2.20)]. Here 6 is 
the angle between the local magnetic field B(r,f) and the di- 
rection to the observer In the case of a random magnetic field 
it is conveni ent to use the local spectral dens ities of the Stokes 
parameters (IGinzburg & Syrovatskiil [T965h . A synchrotron 
map of a source can be constructed as a function of the po- 
sition R± in the plane perpendicular to the observer direction 
from the photons of a fixed frequency. Because of the additive 
property of the Stokes parameters for incoherent photons, we 
can integrate + p*^^ over the line of sight in the observer 
direction weighted with the distribution function of emitting 
particles N{r,E,t) to get intensity I{R±,t,iy). To collect the 
photons reaching the source surface at the same moment, f, 
we performed the integration over the source depth using the 
retarded time t' = t - |r-Rj^|/c as arguments in B(r,f') and 
N{r,E,t'). The integration grid had cell size smaller than 
Lmin- We produce the maps at different frequencies ly and 
for different times to study spectra and variability. Below we 
present synchrotron images simulated with a steady model 
distribution of electrons accelerated by a plane shock of ve- 
locity 2,000 km s"' propagated in a fully ionized plasma of 
number density 0.1 cm"-'. The kinetic model used to simu- 
late the electron distribution was described in detail by ?. The 
magnetic field in the far-upstream region was fixed at 10 ijG 
and it was assumed that magnetic field amplification produced 
a random field of ^ (B^) = lO"'* G in the shock vicinity. The 
effective shock compression ratio is about 4.7. 

The spatial distribution of ^ TeV electrons is highly 
non-uniform and concentrated near the shock due to 
synch rotron losses (see, also IZirakashvili & AharonianI 
l2007l) . The simulated electron energy spectrum at the 
shock front position has asymptotic behavior N (x 7 x 
10-'-[£/TeV]-' '^exp[-(£/7TeV)2] cm-^^TeV"' in the cut- 
off energy regime. 

We note that, in reality, the distribution of the emitting elec- 
trons would be a random function of the position and particle 
energy because of the stochastic nature of both the electro- 
magnetic fields and the particle dynamics. However, such a 
stochastic distribution can currently only be determined with 
particle-in-cell (PIC) simulations and then only for very lim- 
ited energy ranges. The PIC simulation of TeV distributions 
in SNR shocks is beyo nd current computer capabilities (see, 
IVladimirov et al.ll2008l) . 

Our model highlights the effects of a stochastic magnetic 
field on the synchrotron mapping and demonstrates the im- 
portance of the high-energy cutoff region of the electron dis- 
tribution where the synchrotron emissivity depends strongly 
on the local magnetic field. 

4. MAPPING OF A SYNCHROTRON EMISSION SOURCE 

In this section we calculate maps of a synchrotron emission 
source with a random magnetic field. We use a set of parame- 
ters for a stochastic magnetic field sample with y^JW) = lO""* 
G and 5=1.0. Note that a 7 TeV electron has a characteristic 
synchrotron emission energy lii^max ~ 0.1 keV (calculated for 
a uniform magnetic field B = IQ~* G). However, this energy is 
not a good spectral characteristic in a random magnetic field 
because of the presence of the long tails in the probability 



Fig. 2. — Details of Fig[T]in the vicinity of tlie briglit dot Dl. Tlie four panels show, from left to right, ■ I(R^,t, y) maps of synchrotron emission at 0.5, 5, 
20 and 50 keV, respectively. 



distribution function (PDF) of the random field, as is clear in 
Figure [3] In Figure [T] we show ■ I(R± ,t^v) maps at some 
particular time. Shock is propagating perpendicular to the line 
of sight direction. The 5 keV image (left panel a) clearly 
demonstrates the presence of structures: filaments, clumps 
and dots produced by the stochastic field topology. Filaments 
are less prominent at 50 keV (panel b), while dots and clumps 
are present. Polarization maps, where the degree of polariza- 
tion is determined in a standard way also show clumpy struc- 
tures. The degree of polarization reaches very high values in 
high energy clumps. We will discuss in detail the polariza- 
tion properties of synchrotron structures and clump statistics 
elsewhere. 

In Figure |2l the synchrotron emission of a small region 
around the dot Dl is shown for 0.5, 5, 20 and 50 keV. There 
are observable differences in the synchrotron maps indicat- 
ing that some features are bright at high energies and much 
less prominent at lower energies and vice versa. The physi- 
cal reason for this difference comes from the fact that in the 
high-energy cut-off region of the electron spectrum, the lo- 
cal synchrotron emissivity depends strongly on the local field 
value. High order statistical moments of the field dominate 
the synchrotron emissivity at high energies, and even a sin- 
gle strong local field maximum can produce a feature (dot 
or clump) that stands out on the map even after integrating 
the local emissivity over the line of sight. In lower energy 
maps, the contribution of a single maximum can be smoothed 
or washed out by contributions from a number of weaker field 
maxima integrated over the line of sight. The high energy 
map is highly intermittent because the synchrotron emissivity 
depends on high order moments of the random magnetic field 
at the cut-off frequency regime. 

4.1. Spectra and temporal evolution of structures 

In Figure [3] we show the simulated spectral density of the 
surface emissivity of a synchrotron source of size D = 0.3 pc 
with the magnetic field model described above. The spec- 
tra for the regions of different surface emissivity apparent in 
Figure [T] are plotted separately in Figure |3] A dim rectangu- 
lar region A and a bright dot D 1 indicated in Figure [T] have 
drastically different spectra above 0.1 keV in the model. At 
the same time the spectrum of Dl is very similar to the total 
emissivity averaged over the whole map, indicating that the 
total emissivity is dominated by the bright dots and clumps. 

The fact that bright dots and clumps contribute so strongly 
to the total emission means that there will be a clear spec- 
tral difference between the case where a source is filled with 
a uniform magnetic field (i.e., one equal to 0.1 mG in our 
simulation), and the case where the field is turbulent with 
^ (Z?2) =0.1 mG, as we simulate in FiglT] The presence of 
magnetic field fluctuations can substantially increase the syn- 
chrotron emissivity of electrons in the exponential tail of the 
distribution and if this effect is neglected, average magnetic 
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Fig. 3. — Specti'al densities of surface emissivity for regions in Fig.[Tlas 
labeled. The solid curve is the average over the entire region shown in Fig^ 
while the dot-dashed curve is the spectrum produced if the region is filled 
with a uniform magnetic field of fixed magnitude 5=10"'* G, but of random 
direction. 

fields may be inferred that are factors of several larger than is 
in fact the case. 

An important aspect of synchrotron emission in turbulent 
fields is that the emission will vary in time even for a static 
electron distribution. In Figure |4] we show the relative in- 
tensity of the bright unresolved dot Dl [of size ^ 10~^D] as 
a function of time measured in units of Z)/vph, where Vph is 
a characteristic speed of the magnetic turbulence, e.g., the 
Alfven speed. Figure |4] shows that a factor of two or more 
change in intensity will occur on times of a year if Z)/vph ~ 
lO'^sforX -rays with energies > 5 keV. We simulate the im- 
age with minimal pixel size 10~~D = 10'^ cm, while Lmi,, was 
6 X lO'** cm. For SNRs at kpc distances the simulated pixel 
size corresponds to the ~ arcsecond range resolution of Chan- 
dra. The minimal variability time scale depends on the size of 
the structure / if Z > Lmin and the dependence on I is nonlin- 
ear since the photon emission is determined by high statisti- 
cal moments of the random field in the cutoff regime. A one 
year variation-time of Dl at 5 keV requires Vph < 10^ cm s~', 
which is a realistic value for the Alfven velocity for SNRs if 
the magnetic field is > 0. 1 mG. Figure |4] shows that higher 
energy X-rays in bright spots can vary on much shorter time 
scales, another consequence of emission from the exponential 
tail of the electron distribution. Some of the large scale fil- 
amentary structures seen in Fig [T] also show a complex time 
variability. 

5. CONCLUSIONS 

We show that prominent, evolving, localized structures 
(e.g., dots, clumps and filaments) can appear in synchrotron 
maps of extended sources with random magnetic fields, even 
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Fig. 4. — Temporal evolution of dot Dl shown in Fig|2]at 50 keV (dot- 
dashed curve), 20 keV (dashed curve), and 5 keV (solid line). The upper 
panel shows the evolution of Dl at 5 keV on a longer time scale. 

if the particle distribution is smooth and steady. The bright 
structures originate as high-energy electrons radiate effi- 
ciently in local enhancements of the magnetic field intensity. 
The peaks in synchrotron emission maps occur due to high- 
order moments of the magnetic field probability distribution 
function (PDF). 

Even if the PDF of projections of the local magnetic field 
is nearly Gaussian, the corresponding PDF of synchrotron 
peaks, simulated for a spatially homogeneous relativistic par- 
ticle distribution, has strong departures from the Gaussian at 
large intensity amplitudes. This is because of the non-linear 



dependence of the synchrotron emissivity on the local mag- 
netic field in the high-energy cut-off regime of the electron 
spectrum. 

Our model makes two basic predictions. One is that the 
overall efficiency of synchrotron radiation from the cut-off 
regime in the electron spectrum can be strongly enhanced in a 
turbulent field with some ^JIW), compared to emission from 

a uniform field, Bq, where Bq = \J (B^). The second is that 
strong variations in the brightness of small structures can oc- 
cur on time scales much shorter than variations in the underly- 
ing particle distribution. The variability time scale is shorter 
for higher energy synchrotron images. Our estimate of the 
time scales of these intensity variations is consistent with the 
rapid time variability seen in some young SNRs by Chandm. 
The strong energy dependence we predict may be important 
for the future missions NuSTAR, Simbol-X and GRI that will 
image SNR shells up to 50 keV. 

Non-thermal emission in many sources of interest such as 
gamma-ray bursts, AGNs and galaxy clusters originates from 
synchrotron radiation of ultra-relativistic particles in turbulent 
magnetic fields. The effect of random fields on synchrotron 
spectra and time evolution should be accounted for in their 
models. Relativistic flows in GRBs, PWNs and AGNs likely 
produce random magnetic fields with large amplitudes that 
should increase the efficiency of synchrotron radiation from 
electrons/positrons in the spectral cut-off regime. Synchrotron 
radiation in highly stochastic fields also results in energy band 
specific variability of the observed emission and its polariza- 
tion. These specific characteristics of the synchrotron images 
and spectra in the cut-off regime can be used to study the tur- 
bulent magnetic fields in astrophysical sources. 



We thank a referee for constructive comments. A. M. B. 
and Yu. A. U. acknowledge support from RBRF grant 06- 
02-16884 and D. C. E. from NASA grants NNH04Zss001N- 
LTSAand 06-ATP06-21. 



REFERENCES 



Bamba, A., Yamazaki, R., Yoshida, T., Terasawa, T., & Koyama, K. 2005, 

ApJ, 621, 793 
— . 2006, Advances in Space Research, 37, 1439 

Butt, Y. M., Porter, T. A., Katz, B., & Waxman, E. 2008, MNRAS, 386, L20 
Candia, J., & Roulet, E. 2004, Journal of Cosmology and Astro-Particle 
Physics, 10, 7 

Casse, F, Lemoine, M., & Pelletier, G. 2002, Phys. Rev D, 65, 023002 
Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 
Ginzburg, V. L., & Syrovatskii, S. I. 1965, ARA&A, 3, 297 
Patnaude, D. J., & Fesen, R. A. 2008, ApJ (submitted) ArXiv 0808.0692 



Pohl, M., Yan, H., & Lazarian, A. 2005, ApJ, 626, LlOl 

Uchiyama, Y, Aharonian, F. A., Tanaka, T., Takahashi, T., & Maeda, Y. 2007, 

Nature, 449, 576 
Vink, J., & Laming, J. M. 2003, ApJ, 584, 758 

Vladimirov. A. E., Bykov, A. M.. & Elhson, D. C. 2008, ApJ, 688, 1084 
Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2005, A&A, 433, 229 
Zirakashvili, V. N., & Aharonian, F 2007, A&A, 465, 695 



